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The eigenstate thermalization hypothesis (ETH) attempts to bridge the gap between quantum 
mechanical and statistical mechanical descriptions of isolated quantum systems. Here, we define 
unbiased measures for how well the ETH works in various regimes, by mapping general interacting 
quantum systems on regular lattices onto a single particle living on a high-dimensional graph. By 
numerically analyzing deviations from ETH behavior in the non-integrable Ising model, we propose 
a quantity that we call the n-weight to democratically characterize the average deviations for all 
operators residing on a given number of sites, irrespective of their spatial structure. It appears 
to have a simple scaling form, that we conjecture to hold true for all non-integrable systems. A 
closely related quantity, that we term the n-distinguishability, tells us how well two states can be 
distinguished if only n-site operators are measured. Along the way, we discover that complicated 
operators on average are worse than simple ones at distinguishing between neighboring eigenstates, 
contrary to the naive intuition created by the usual statements of the ETH that few-body (many- 
body) operators acquire the same (different) expectation values in nearby eigenstates at finite energy 
density. Finally, we sketch heuristic arguments that the ETH originates from the limited ability of 
simple operators to distinguish between quantum states of a system, especially when the states are 
subject to constraints such as roughly fixed energy with respect to a local Hamiltonian. 


I. INTRODUCTION 


Statistical mechanics and quantum mechanics have 
been the cornerstones of modern physics for nearly a cen¬ 
tury. Both formalisms have been put to test in a wide 
variety of scenarios, and both have invariably given de¬ 
scriptions that are accurate as well as consistent with 
each other within their regimes of validity. Their starting 
points, however, are fundamentally different. Quantum 
mechanics describes unitary evolution of isolated systems 
via wave functions or pure states, whereas statistical me¬ 
chanics ascribes microcanonical ensembles - which are 
inherently mixed states - to them. Why, then, do the 
two theories concur? 

A crucial step towards answering this question was 
taken by Berry who conjectured that the eigenstates 
of quantum Hamiltonians whose classical counterparts 
are chaotic behave as if they were drawn randomly from a 
Gaussian distribution. In other words, they resemble ran¬ 
dom superpositions of classical conhgurations at the same 
energy density in the thermodynamic limit. Deutsch 
argued that Berry’s conjecture holds for integrable sys¬ 
tems perturbed away from integrability. The argument 
stems from the intuition that integrability-breaking per¬ 
turbations allow quantum systems to mimic the central 
maxim of classical statistical mechanics. On the one 
hand, classical ergodic systems explore all accessible mi¬ 
crostates with equal likelihood within the time scale of 
typical measurements and hence, justify equating time 
and ensemble averages. On the other, quantum systems 
can access all classical “microstates” at once from a single 
eigenstate, so no time or ensemble-averaging is required 
and each eigenstate effectively resembles a microcanoni¬ 
cal ensemble. Srednicki Q then explicitly derived statis¬ 
tical distributions for quantum particles assuming only 
Berry’s conjecture, and postulated the eigenstate ther¬ 
malization hypothesis (ETH). The ETH states. 


In ergodic quantum systems, eigenstates at 
finite energy density give rise to expecta¬ 
tion values for “simple operators” that vary 
smoothly with energy, with fluctuations that 
are exponentially suppressed in the system 
size. The off-diagonal matrix elements of 
these operators between nearby eigenstates 
also vanish exponentially. 

Here, “simple operators” usually refers to operators that 
involve a very small number of degrees of freedom com¬ 
pared to the system size, although the precise definition 
is given very rarely [^. 

The ETH is only believed to hold for simple operators 
and states with hnite energy density in the thermody¬ 
namic limit. With increasing complexity of operators, 
their ability to tell eigenstates with nearby energy apart 
improves. The extreme case is that of a projection opera¬ 
tor pn = |n.)(n| onto a given eigenstate |n); this typically 
has support on all sites of the system and it obviously dis¬ 
tinguishes \n) from all other eigenstates perfectly. There¬ 
fore, a refined statement of the ETH demands the defini¬ 
tion of quantitative measures for how distinguishable two 
states are when a restricted set of operators is measured. 
Moreover, we will see that choosing this restricted set to 
be all operators of certain complexity (rather than some 
arbitrarily chosen operators) permits the definition of a 
“basis-independent measure” to systematically quantify 
how well the ETH works and how distinguishable nearby 
states are as a function of the system size, eigenstate en¬ 
ergy and operator complexity. 

In this work, we define two such measures: 


1. the n-weight of a Hermitian operator A, 
P„(A) = ^|Tr(AOl"))|' 

i 


( 1 ) 
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where for all t form a complete orthonormal 
basis (with respect to the Hilbert-Schmidt norm) 
of operators that have support on n sites, and 

2. the n-distinguishability of two operators Ai and 

^2, 


0n{Ai,A2) = cos 


_1 ^rt(^l) + Pn{A2) — Pn{Ai — A 2 ) 


2^P„iA,)P„iA2) 


( 2 ) 


These measures are defined for a lattice system with the 
Hilbert space IH = ]([® Hi, Hi being the Hilbert space of 
the site. The intuition behind these quantities is an 
alternative view of quantum many-body dynamics. As 
will be described in Sec. im Heisenberg evolution of a 
many-body operator can be mapped to a quantum me¬ 
chanical problem of a single particle hopping on a high¬ 
dimensional graph. Each site of the graph represents a 
basis operator, e.g., a direct product of Pauli operators 
if each lattice site hosts a single qubit. The "particle 
wavefunction" is determined by the overlap of operator 
A with the basis operator. The graph has a layered struc- 
ture with all n-site operators O) ' in the n-th layer. In 
this picture, the n-weight Pn{A) is simply proportional 
to the probability of the particle being in the n-th layer 
of the graph. The n-distinguishability of two operators 
Ai and A 2 is the angle between the projections of the 
corresponding vectors in the Hilbert space of operators 
onto the subspace of the graph sites that constitute the 
n*^ layer. A large angle means Ai and A 2 are easier 
to distinguish with n-site operators. 

Most previous studies have numerically analyzed the 
expectation values or off-diagonal matrix elements of a 
select set of simple operators dwell ing in contiguous re¬ 
gions of space either at long times pulll or in different 
eigenstates [iM3. Infinite temperature studies are the 
only exceptions that we are aware of that include non¬ 
contiguous regions of space (Tsl - f^ . In contrast, our mea¬ 
sures treat all operators with support on a given number 
of sites on equal footing, irrespective of whether the sites 
are adjacent or spatially separated. If A is chosen to 
be a density operator, then the n-weight turns out to 
be closely related to the second Renyi entropy of regions 
with size n and smaller. However, the Renyi entropy of 
an n-site region has contributions from fewer site oper¬ 
ators too, while the n-weight is a direct measure of the 
amplitude of n-site correlation functions in this density 
operator. 

Using the two measures we define, we study two re¬ 
lated quantities that characterize the difference between 
two neighboring eigenstates, with density operators pi 
and p 2 - The first quantity is the n-weight of the differ¬ 
ence, Pn {pi — Pi)- By performing numerical exact di- 
agonalization on a simple non-integrable model, an Ising 
model with parallel and transverse fields, we discover that 
Pn {pi — P 2 ) attains a simple scaling form as a function 
of n and the density of states in the part of the spectrum 
from which the two states are drawn. We conjecture 


that a similar form holds in all non-integrable systems. 
Our result on (pi — P 2 ) also clarifies a rather pecu¬ 
liar feature of quantum ergodic systems. According to 
the ETH, one would naively think that a complicated 
(high n) operator has a better chance of distinguishing 
two neighboring eigenstates than a simple operator. We 
discover that surprisingly, the opposite is true for typical 
operators: simple operators on average are better at dis¬ 
tinguishing between nearby eigenstates than complicated 
operators are. In fact, there is a critical size of operators, 
which only depends on the system size and the Hilbert 
space dimension at each site, at which measuring a ran¬ 
dom operator reveals no information about the system. 
This critical operator size occurs well beyond half of the 
system size, and was consequently missed by previous 
studies that focused on simple operators. 

The second quantity we study is the n- 
distinguishability for the same pair of eigenstates, 
dn{pi, P 2 ), which helps to reconcile the intuition and 
the numerical observation mentioned above. This angle 
is found to be small when the vectors are projected 
onto a low-n subspace, determined by simple operators, 
but to saturate to 7r/2 for more complex operators. 
Both these behaviors are found to be strikingly different 
in integrable systems. The large n-distinguishability 
suggests that although most large operators are poor 
at distinguishing between neighboring eigenstates, there 
are a few fine-tuned ones that are proficient at it. Thus, 
the optimal n-site operator can differentiate between the 
two states better if n is large, but a randomly chosen one 
works better for smaller n. The combination of the two 
measures and thus provides a more systematic and 
accurate understanding to the nature of thermalization. 

The remainder of the paper is organized as following. 
In Sec. |H]we describe the Fock space of operators and set 
the general framework of our approach. The two quan¬ 
tities n-weight and n-distinguishability are defined and 
their general properties are discussed in Sec. IHII Sec. 
m is devoted to analyzing numerical results in the Ising 
model, while Sec. El is dedicated to further discussions 
and conclusion. 


II. FOCK SPACE OF OPERATORS 


We start by defining a mapping from a lattice many- 
body system to a single-particle quantum mechanics 
problem on a high-dimensional graph. For concrete¬ 
ness and simplicity, let us choose the many body sys¬ 
tem to be a spin-1/2 model with Hamiltonian H on a, 
chain with L sites. For a spin-1/2 model, the Hilbert 
space of site r is spanned by an orthonormal basis con¬ 
sisting of the identity matrix and the Pauli matrices: 

(Ir,cr/, ct/,ct/) = = 0,1,2,3, normalized as 

Ti (^0(^r,i)0(^r,j)) = 5ij. They satisfy the orthogonality 
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condition 

3 1 . 

E (^lr,^)Olr,i) = SadSbc-^6abS,d = = X, - ^ (3) 

i=l 

where D = 2 is the dimension of the Hilbert space at 
each site, = SadSbc is the “swap” operator acting 

on site r, whose expectation value computed over two 
copies of any state p is related to the second Renyi en¬ 
tropy of site r in that state [ 1 ^: Sr = — logTr^ (Tr^-p)^ = 
— logTr [(p 0 p)Xr], and Wr is its traceless part. Here, 
r denotes all sites except site r. The Hilbert space 
of the whole system is spanned by multisite operators 

^{rl,,ic\a=l...n} ~ ^SiM) ® ^(*'2 ,*2) ■ ' • ,in) = \ 

where n represents the number of sites on which 
is one of the Pauli matrices, and sites absent in the set 
{ra} are assumed to host the normalized identity op¬ 
erator in Henceforth, n will be referred to 

as the size of the operator The subscript i in¬ 

dexes all operators with size n, so it runs from 1 through 

{D"^ — 1)" ^ ^ ^ We will refer to operators 

with small (large) n as simple (complicated) operators. 

Now, we construct a graph by assigning a node to each 
basis operator 0\ . The nodes are sorted by operator- 
size, so the layer of the graph is comprised of nodes 
corresponding to n-site operators, as depicted in Fig. [T] 
The total number of operators, or nodes in the graph, is 
, while the number of nodes in each layer is fnD^^. 
The outline in Fig. [T] depicts the number of nodes in 
each layer, fnD^^. For L ^ D^, fn is maximized at 
n = n* = (1 — 1/D^)L which equals 3L/4 for D — 2. It 
is easy to check that n* is the mean, median as well as 
the mode of /„. Physically, this means that a randomly 
chosen basis operator is most likely to have a size of n*, 
since each site has — 1 non-trivial operators and a 
single trivial one. 

The orthogonality of the operator basis allows us to ex¬ 
pand an arbitrary many-body Hermitian operator A{t) in 
this basis as A(t) = ^ with real coeffi¬ 

cients: 

V;„r(A;t) =Tr(24(t)0f)) (4) 

The vector %l!ni{A\t) can be viewed as a "single-particle 
wavefunction" of a particle hopping on the graph we de¬ 
fined. There are two key requirements for interpreting 
ipniiA', t) as a sensible single-particle wavefunction within 
first quantization. Firstly, the total probability density 
of the particle must be conserved; this is guaranteed be¬ 
cause I form an orthonormal basis for operators in 

the Hilbert space, which implies 

ElV'nr(R;t)|"=Tr(H 2 ) ( 5 ) 



Figure 1: Schematic of the tree network, constructed by ar¬ 
ranging operators according to their size. The identity opera¬ 
tor is on the extreme left, followed by the single-site operators, 
and so on. The total number of operators is , while the 

fraction in the layer is /„ = ^ n ) ’ max¬ 

imum for n = (1 — 1/D^)L — |L = n* for D = 2. Note that 
fn grows exponentially with n for small n, so the apparent 
linear growth in the number of dots with n in the figure is for 
ease of depiction and should not be taken literally. 


a manifestly invariant quantity under unitary time- 
evolution of A. Secondly, it must satisfy Schrodinger’s 
equation. Indeed, A{t) follows Heisenberg time evolu¬ 
tion: A{t) = i[H,A{t)], so the time-evolution of ipne is 
given by 


= Ar([A{t),H]0f) 

= R?/V’nT'(R;t) (6) 

where the hopping matrix element between two nodes 
and ^ is defined as 


= Tr ( 


-H, O 




_ njn n _ { nj 


[m?') 


( 7 ) 

Thus, ifniiA'A) is a sensible wavefunction, and its dy¬ 
namics are governed by the hopping Hamiltonian 0. 
Crucially, if the physical Hamiltonian contains at most 
fc-spin terms, Wlf) vanishes for |n — n'| > k and there¬ 
fore, satisfies a notion of locality on the graph. Thus, 
the time-evolution of a general operator in a many body 
system has been recast into the problem of a single par¬ 
ticle governed by a local hopping Hamiltonian on a high¬ 
dimensional graph. The “energies” of the particle con¬ 
sists of all possible differences Ei — \ 

between pairs of energies of H, so the single-particle spec¬ 
trum that results from a generic El with no degeneracies 
is particle-hole symmetric with zero eigenvalues and 
7 ^ 2 L _ jjL jion-zero ones. The above construction is rem¬ 
iniscent of mappings of states in Fock space to a Cayley 
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tree, which allows one to view integrable systems with 
local conservation laws in real [3 momentum 
space as a localized particle on a suitably defined tree. 
Our construction, on the other hand, describes the Fock 
space of operators instead of that of states, and thus is 
closer to the approach adopted by Ref. (2^1 for construct¬ 
ing integrals of motion to the describe the many-body 
localized phase. 

The graph construction facilitates extracting different 
kinds of information about quantum ergodicity with dif¬ 
ferent choices of A. If we choose A to be the density 
matrix p for an eigenstate that satisfies the ETH, then 
'4’neiA) is expected to coincide with its thermal value, 
determined only by the energy density in that state, for 
small n, but can be depend on other details for larger 
n. On the other hand, picking A to be the difference 
between two density matrices tells us how various corre¬ 
lators differ in the two states, and what kind of observ¬ 
ables must be measured in order to distinguish between 
them. Finally, if A is a physical observable O, its time- 
evolution tells us how this physical quantity evolves into 
a superposition of other operators with time. In partic¬ 
ular, suppose we start from a small-n operator and let it 
evolve in time. In general, it will evolve into a superpo¬ 
sition containing many large-n operators. Equivalently, 
one can think of (O; t) as an infinite temperature cor¬ 
relation function, so its time evolution tells us how small 
operators develop correlations with large operators over 
time. In the language of the hopping particle, this cor¬ 
responds to the particle starting near the low-n end of 
the graph and spreading towards larger n nodes. Thus, 
chaotic behavior of operators turns into delocalization of 
the particle on the graph. 


III. n-WEIGHT AND n-DISTINGUISHABILITY 

A central quantity that we will work with in this paper 
is the n-weight of an operator A, Pn{A), defined as 

P^{A)=J2\MA)\^= E 

( 8 ) 

where the sum over I runs over all possible choices of 
n-site operators. In the graph picture of Sec. m this 
is the total probability density in the layer of the 
graph. Physically, P„(A) tells us how complex the cor¬ 
relators one must measure in order to reconstruct A are. 
Computing it directly is computationally taxing, as it 
entails computing traces, one for each operator in 
the Hilbert space. Fortunately, the computation can be 
simplified via a generating function, as follows. 

Performing the sums over {ia} for fixed {tq,} and using 
Eq. ([3]) gives 


^ E n Tr [(^ ® (9) 

R-n r£Rn 


where i?„ is a region of size n (not necessarily connected), 
and the sum is over all n-site regions. Now we define 
a generating function 


FiA;z) = J2^^Pn{A) 


( 10 ) 


n =0 

which can be explicitly written as 

L 


F{A;z) = ^Tr 


A O A P (1,. -I- DzWr 


( 11 ) 




Here the sum is over all regions R, composed of kR 
sites and R denotes the complement of R. Pn{A) is then 
determined by Fourier transforming P(A; e*"®): 


Pn{A) = 


27r 


/' 


d6»e-“^P(A;e™®) 


( 12 ) 


Since n £ [0, A] is linear in system size, it is actually suf¬ 
ficient to calculate F(A; e™®) for the L + 1 discrete values 
of 0 w = 0,1,.., L. Therefore we have trans¬ 

lated the calculation of Pn{A) for all A to L -I- 1 operator 
trace computations in the doubled Hilbert space. Alter¬ 
natively, we can also use the second line of Eq. m and 
calculate Pn by performing partial trace calculations 
on a single copy of the system. 

When A is a density matrix p, the generating function 
also makes it transparent that (A) is related to the sec¬ 
ond Renyi entropy of a region R, Sr = — logTr (pr) = 

-logTr From®, we have 


Pn{p) = 


1 

pL-n 


E E 

Rn RkQ^Rn 





(13) 


The sum over R„ and Rk can be combined into a single 
sum by introducing suitable combinatorial factors. Defin¬ 
ing e~^’‘ = i e., the average of e~^^k over all 

A:-site regions, we get 



Thus, there is a simple relationship between P„, the 
“single particle density on the graph”, and entanglement 
properties of the many body state. 

Based on the n-weight defined for each operator A, we 
define a second quantity, the n-distinguishability between 
two operators Ai and A 2 , as 


6 >„(Ai, A 2 ) = cos 


-1 Pn(Ai) + P„(A2) — P„(Ai — A 2 ) 


( 15 ) 
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Pi 



Figure 2: Schematic illustration of the behavior of neighbor¬ 
ing eigenstates pi and p 2 on being projected onto the space 
of small-n (left) and large-n (right) operators, p^"^ is short¬ 
hand for the projections tpnipi) defined in the text. The three 
directions together depict the full Hilbert space of operators, 
while the horizontal plane represents its projection onto the 
space of operators of size n. pi and p 2 are mutually orthog¬ 
onal vectors in the full space. However, they appear nearly 
identical when projected onto simple operators as shown on 
the left, but look quite different for larger operators as shown 
on the right. The angle 6n will be calculated numerically for 
the non-integrable Ising model in Sec. m 


On (Ai,A 2 ) is simply the angle between the two vectors 
ifn{Ai) = tpne{Ai) = and tfn{A2) defined 

similarly, i.e. the two “single particle wavefunctions" cor¬ 
responding to Ai and A 2 , projected to the graph sites 
corresponding to size-n operators. The angle On thus 
measures how different the two operators are if only n-site 
operators are measured. A small implies that Ai and 
A 2 look similar in all size-n measurements, while a large 
On ^ means Ai and A 2 can be easily distinguished by 
n-site operators. This is sketched in Fig. [21 where we 
have chosen A 1.2 to be two density matrices pi 2 in an¬ 
ticipation of the discussion in the next section. In Sec. 
IIV Bl we will use On to distinguish between neighboring 
eigenstates and show that indeed, they appear similar for 
simple observables and different for complicated ones. 


IV. EIGENSTATE THERMALIZATION IN THE 
ISING MODEL 

In this section we apply the new measures we define 
to study eigenstate thermalization in a prototypical non- 
integrable spin model, namely, the ID Ising model with 
transverse and longitudinal fields, given by 

( 16 ) 

r 

H is integrable if any one of J, and hz vanishes, 
but is non-integrable otherwise. We choose J = 0.5, 
hx = —0.74 and hz = 0.35 as the non-integrable parame¬ 
ters, and J = 0.5, hx = 0.35 and = 0 as the integrable 
ones. Open boundary conditions and the extra term on 
the first site, hzcrf, ensure that translation and inversion 



Figure 3: Pn{Ap) vs n (above) and Pn.{Ap)/fn vs n (below) 
for ~ 200 randomly chosen pairs of neighboring eigenstates for 
L = 14 sites for the ergodic (left) and the integrable (right) 
Ising model. Here, Ap — pi — p2 is the difference between the 
density matrices of the two eigenstates. The color is propor¬ 
tional to the density of states, with blue (red) representing 
states in regions of the spectrum with low (high) density of 
states. The black dashed line marks Pn{Ap) — 2/„. 


symmetries are broken so that there are no conserved 
quantities in the non-integrable case. This is unlike sev¬ 
eral recent works which retained translational symm etry 
and hence, conserved the total momentum [ 3,0 [ 13 ; l27l |. 
The energies and eigenstates are obtained by exact diag- 
onalization of systems of upto L = 14 sites. 


A. Gomparison of eigenstate n-weights 

As stated in the introduction, the ETH says that the 
expectation values of simple operators are equal in nearby 
eigenstates of chaotic Hamiltonians, up to exponentially 
small corrections in the system size. This automatically 
ensures that each eigenstate resembles a “microcanonical 
ensemble”, i.e., an equal admixture of nearby eigenstates, 
in the thermodynamic limit and hence yields the ETH as 
stated in the introduction. Thus, we first compare pairs 
of neighboring eigenstates pi and p 2 by computing the 
total squared difference in the expectation values of all 
operators of size n, 

P„(Ap) = ^((0S">)^_-(0<">)j' (17) 

where Ap = pi — P 2 , and study its dependence on n 
and the energy of the pair. As shown in the upper pan¬ 
els of Fig. 131 this quantity has the anticipated behav¬ 
ior for small n: it increases with n and is larger when 
the density of states is lower. A closer inspection, how¬ 
ever, reveals that the n-dependence seen here is decep¬ 
tive, and cannot be used to declare eigenstate thermal- 
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ization. In particular, the curves approximately follow 

the fraction of operators of size n, fn = 

upto an overall proportionality constant. In fact, the 
“infinite temperature” eigenstates - states near the mid¬ 
dle of the spectrum where the density of states is high¬ 
est - have Pn{Ap) fa 2/„ = Tr [(Ap)^] /„. Moreover, 
the n-dependence is roughly the same even for the in- 
tegrable Ising model. Thus, we conclude that the bare 
n-dependence of P„(Ap) is primarily determined by the 
number of operators of size n, not by the integrability 
properties of the Hamiltonian. 

Therefore we study the average density per site 
Pn{Ap)/fn - the mean squared difference in the expec¬ 
tation values of size-n operators between neighboring 
eigenstates, upto an overall proportionality constant of 
. As is shown in the lower panels of Fig. |3l the 
n-dependence of this quantity is clearly different for er- 
godic and integrable systems. However, its n-dependence 
for the ergodic system is the exact opposite of what one 
would naively expect from ETH. Indeed, the lower panels 
of Fig. [3] show that on average, large operators are actu¬ 
ally worse at distinguishing between neighboring eigen¬ 
states than small operators are, irrespective of whether 
the Hamiltonian is integrable or not. 

The fact that P„(Ap)//„ decreases with n for small n 
simply says that on average, simple operators store more 
information about the state of the system compared to 
complicated ones. For integrable systems as well as for 
ground states of ergodic systems, this statement is eas¬ 
ily understood because nearby eigenstates can be distin¬ 
guished by simple operators. Fig. |3] says that random 
simple operators can split finite energy density states 
of ergodic Hamiltonians as well, but the efficiency with 
which they can do so decreases with increasing density 
of states. For the infinite temperature states (i.e., states 
at the part of the spectrum with the largest density of 
states) Pn/fn is almost independent of n, which means 
that there is no difference between simple and compli¬ 
cated operators, since the state is essentially a random 
state in the Hilbert space. 

A peculiar feature of Fig. [3] is that random operators 
of size n = n* fa 3L/4 = (I — 1/D‘^)L cannot distinguish 
between any pair of states, irrespective of the Hamilto¬ 
nian. This can be understood heuristically as follows. 
Since each site has — 1 non-trivial operator and a sin¬ 
gle trivial operator, a random operator has size n*. The 
concurrence of Pn{Ap)/fn curves a.t n = n* reflects the 
fact that measuring a random operator does not reveal 
any information about the state of the system. Random 
operators with n > n* can again distinguish between 
neighboring eigenstates. Unlike simple operators, how¬ 
ever, the efficiency with which they can do so in ergodic 
systems increases with increasing density of states. In 
other words, random operators with n < n* are better 
at splitting low lying excitations than at splitting finite 
energy density states, whereas operators with n > n* are 
better at the opposite. 



Although Pn/fn decreases with n for both ergodic 
and integrable systems, its dependence on the density 
of states is clearly different, as can be seen from Fig. |31 
Thus, we highlight the difference by plotting the same 
data as a function of the density of states g{E) in Fig. U) 
Interestingly, for the non-integrable model we find that 
both the n-dependence and the density-of-states depen¬ 
dence can be captured by a clear scaling form given by 

(18) 

In \ 9 oo ) \Ij/ 

where goo = maxgi(E), a is a positive constant and 
S{n/L) is a scaling function of 0(1) that depends only 
weakly on n/L] for the data shown in this section it is 
the exponential of a simple polynomial (See Fig. 01 right 
bottom). Eq. (fT31) depends mainly on generic properties 
of the system such as the density of states of the spec¬ 
trum, the Hilbert space dimension on each site and the 
system size. These properties are obviously common to 
both integrable and ergodic systems; in fact, the scaling 
of g{E) with system size shows no distinction between 
them, as one can see in the left panel of Fig. 01 However, 
we see that integrable systems do not have simple scaling 
form as (ITSl) for P„. We hypothesize that the scaling be¬ 
havior in Eq. (|18ll is a generic property of non-integrable 
systems, and hope that this conjecture can be tested in 
other systems. 


B. n-distinguishability of neighboring eigenstates 

How do we reconcile the decrease of Pn{Ap)/fn with 
n in Fig. 131 with the anticipation from ETH that ex¬ 
pectation values of large operators, in some sense, de¬ 
viate more than those of small operators? In order to 
resolve this counterintuitive behavior, we compute the 
n-distinguishability 9nipi, P 2 ) for the same pair of neigh¬ 
boring eigenstates pi and p 2 and present the result in Fig. 
[5l As we discussed earlier, is the angle between two 
vectors in the size-n Hilbert space ifn{pi) and ifn{p 2 ), 
and the vectors are lists of average values of all size-n 
operators in the two states pi , p 2 respectively. 

As is shown in Fig. 9n increases monotonically 

with n for most pairs of states in the ergodic system 
and quickly saturates at the maximal value 7r/2 before 
n reaches L/2. By combining this observation with the 
behavior of Pn/fn observed earlier, one can understand 
better what happens with increasing n. (pi — P 2 )/fn = 
\'f’n{pl)/^/f/i.-i’n{p 2 )/^/]n\'^ is the 2-norm squared of the 
difference between the two vectors, i.e., the Euclidean dis¬ 
tance between them. Although this distance decreases 
with increasing n, the decrease is mainly due to the 
shrinking of the norm of each vector, and the angle be¬ 
tween them is actually increasing. For large n the two 
vectors are both very short (which means a typical size-n 
operator has a small average value), but they are almost 
always exactly perpendicular. In contrast, at small n the 
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I Density of stat^ 


Pn{Ap) / fn vs density of states 


Scaling of non-integrable fits 



Figure 4: Left: j- log [gaa/g{E)] vs the energy density EjL for various system sizes for the ergodic (above) and the integrable 
(below) Ising model. Except near the band edges, the curves are indistinguishable, indicating that g{E)/goo grows exponentially 
with L with an energy density dependent exponent. Middle: log [Pn{Ap) / f„] vs log [gcx/g{E)] for various n for L = 14 for the 
ergodic (above) and the integrable (below) Ising model, and straight line fits to the ergodic data. The fits are good for n > 3. 
A similar fitting procedure for other system sizes yields the panels on the right, where we show that the slopes (above) and the 
intercepts (below) of the straight lines are simple functions of n/L. 


angle is small, meaning the simple-operator average val¬ 
ues in the two neighboring eigenstates are well-correlated. 
In other words, as n increases, although a randomly cho¬ 
sen size-n operator does a worse job distinguishing the 
two states, there exists a particular choice of operator 
which can distinguish the two states better. When the 
two vectors become perpendicular at large n, the two 
states can be distinguished completely if we simply use 
the n-size operator (Ap), i.e., the projection 

of Ap to the n-size subspace. 

At the maximum of the density of states, a.t E/L — 0, 
we see a peak of for small n. The ETH is expected to 
work best in this regime, but the behavior of indicates 
that it is violated dramatically. However, this feature 
can be safely ignored because at infinite temperature, the 
simple operators’ average values are almost vanishing so 
the angle between them is inconsequential. 


V. CONCLUSION AND DISCUSSIONS 

In summary, we proposed two related new measures 
of thermalization, which help uncover a finer structure of 
the process of thermalization than was previously known. 
By mapping the Heisenberg evolution of operators in a 
quantum many-body system onto the Schrodinger evo¬ 
lution of a single particle on a high dimensional graph. 


the complexity of an operator can be measured by its 
probability distribution in the space of all operators. We 
defined the n-weight of an operator, which measures its 
weight in size-n operators. We studied the behavior of 
the n-weight for the difference between the density op¬ 
erator of two neighboring eigenstates in the energy spec¬ 
trum, denoted as P„(Ap), which tells us how different 
these two states are if we only measure size-n operators. 
We found that Pn{Ap) in the non-integrable Ising model 
follows a simple scaling behavior as a function of the size 
n and the density of states. In contrast to naive ex¬ 
pectation, we saw that large operators on average did 
a worse job at differentiating between two neighboring 
eigenstates in a non-integrable system. In particular, 
there is a critical operator size, which is 3/4 of the sys¬ 
tem size for spin-1/2 models, at which operators on aver¬ 
age cannot distinguish between any states because they 
themselves are “random operators”. This counterintuitive 
conclusion was explained when we investigated the other 
measure of thermalization, the n-distinguishability, de¬ 
fined as an angle between two Hilbert space vectors, each 
corresponding to the projection of an operator onto the 
subspace of size-n operators. We saw that the angle was 
small for low n and saturated to 7r/2 for large n, which 
suggests that two operators (which are orthogonal in the 
whole space) are indeed almost orthogonal in the sub¬ 
space of size-n with large n, while they are almost paral- 
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Figure 5: The angle between the projections of two neigh¬ 
boring eigenstates onto the space of operators of size n, On 
as a function of energy density for various n for L = 14 for 
the ergodic (above) and the integrable (below) Ising model. 
For small n, On is small for most states in the ergodic model, 
but is large in the integrable one. For large n, it is very close 
to 7r/2 in both figures, thus proving that neighboring eigen¬ 
states are nearly orthogonal when projected onto the space 
of n-site operators with large n. Large On at small n and 
large g{E) {g{E) is large near the middle of the spectrum and 
small near the edges; see Fig. 0 is due to the vectors 
and becoming very small in magnitude; then the angle 
between them is meaningless. The sharp change in On at the 
lowest and the highest energies in the ergodic case originates 
from the fact that these states do not satisfy the ETH at the 
current system size. 


lei for small n. In other words, this means that although 
a generic size-n operator does a poor job distinguishing 
two nearby eigenstates for large n, there exists a par¬ 
ticular choice of size-n operators which can distinguish 
them almost perfectly. For small size operators, a ran¬ 
dom operator has an average value that’s more different 
in these two states than a random large size operator, 
simply because the expectation values increase with de¬ 


creasing n. However no small operator can completely 
distinguish the two density operators since the two cor¬ 
responding vectors are almost parallel. In short, random 
operators are better distinguishers for smaller n, but the 
best distinguishers are large operators. 

The “particle-on-the-graph” picture also allows us to 
understand, heuristically, why the ETH holds in the first 
place. It follows from the fact that for vectors living 
in a high-dimensional space (which in this case are the 
states of the particle in its Hilbert space), measuring just 
a few components of two vectors is insufficient to de¬ 
duce whether they are orthogonal or not, especially if 
the components along some of the measured directions 
are almost equal. The role of the Hamiltonian is twofold: 
(i) being ergodic, it ensures that the entire space of op¬ 
erators is connected, as opposed to integrable systems 
where it effectively decouples into subspaces with fixed 
values of the conserved quantities (ii) being a sum of sim¬ 
ple operators, it approximately restricts the expectation 
values other simple operators can take in nearby eigen¬ 
states, i.e., eigenstates with the same energy density in 
the thermodynamic limit. 

The intuition about the role of Hamiltonian can be fur¬ 
ther clarified as follows. A density operator p = \4>){(l>\ of 
an eigenstate can be determined by the following equa¬ 
tions: 

[p,H] = 0 

Ti{pH’^) = E'^, m=l...D^-I (19) 

All density operators for eigenstates can be written as 
a linear superposition of the identity operator and H, 

. If one has access to all the coefficients of 
this expansion, then one can tell that two different eigen¬ 
states are mutually orthogonal. However, the coefficients 
of the low powers of H are slowly varying functions of 
E, so are almost identical for eigenstates with nearby 
energies. Therefore, if we restrict the measurements to 
operators captured by small powers of iL, then it is clear 
why the nearby eigenstates appear almost identical. The 
more nontrivial statement that requires non-integrability 
of the Hamiltonian is the fact that the expectation values 
of other simple operators also depend on the energy of an 
eigenstate smoothly, probably because the high powers of 
H have a negligible contribution to simple operator av¬ 
erage values. More rigorous analytic and numerical work 
is required to make this discussion more complete, which 
we will leave for the future. 

Phase space based reasons for why isolated quantum 
systems fulfil the basic tenets of statistical mechanics 
have been suggested in the past to argue that most sys¬ 
tems, at most times and for most observables, behave 
as if they belonged to a thermal ensemble (28l - [3(]j . The 
proofs there were based on Levy’s lemma which is the 
analog of the law of large numbers for vectors in a high¬ 
dimensional space. 

All the results in this work concern expectation val¬ 
ues of operators in eigenstates. The off-diagonal ma¬ 
trix elements of an operator between different eigenstates 
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{n\0\m) are responsible for time-evolution properties of 
the system. Whether thermalization occurs in time evo¬ 
lution of a state that is not an eigenstate, and how fast 
it occurs, depends on two properties, the off-diagonal 
matrix elements of observables, and the energy differ¬ 
ence between eigenstates. As long as the energy dif¬ 
ferences are sufficiently incommensurate, time evolution 
causes decoherence and thus leads to equilibration (2^ 
Is^ . However, different operators may thermalize at dif¬ 


ferent rates, which is an interesting problem that can be 
studied using the new measures we defined. 
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